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O Abstract 

We review theory of the leading twist nuclear shadowing, and describe phe- 
nomenon of post-selection suppression of leading parton spectrum (effective 
fractional energy losses) in the proximity of the black disk regime. We argue 
that 2 — > 2 mechanism dominates in the inclusive leading pion production in 
d-Au collisions and explain that the post-selection naturally explains both 
the magnitude of the suppression of the forward pion production in d-Au 
collisions and the pattern of the forward - central correlations. At the same 
time this pattern of correlations rules out 2 — >■ 1 mechanism as the main 
source of the inclusive leading pion yield. It is demonstrated that the mech- 
anism of the double parton interactions gives an important contribution to 
^ the production of two leading pions in pp scattering opening a new way to 

study correlations of leading quarks in the nucleon. The same mechanism is 
enhanced in dAu — > n°n + X collisions and explains the dominance of Aip 
independent component and suppression of the away side peak. 

Keywords: leading twist shadowing, black disk regime, forward hadron 
production 
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1. Theory of the leading twist nuclear shadowing 

The connection between nuclear shadowing and diffraction was estab- 
lished a long time ago pQ. One can understand it as a manifestation of the 
unitarity as reflected in the Abramovsky, Gribov, Kancheli (AGK) cutting 
rules [2]. Accuracy of the theory for the hadron - nucleus interactions is 
on the level of few % which reflects small admixture of non-nucleonic de- 
grees of freedom in nuclei and small off-shellness of the nucleons in nuclei 
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as compared to the soft strong interaction scale. The Gribov logic was suc- 
cessfully applied for the description of the 7 A total cross sections. However 
an additional step is necessary to calculate the nuclear parton distributions 
at small x [3J. One has to combine unitarity relations for different cuts of 
the diagrams corresponding to the different final states for the interaction 
of the hard probe with a nucleus, with the QCD factorization theorem for 
hard diffraction |4j. In the limit of thin target when only double scattering 
is important one obtains 

A [xf j/A (x, Q 2 , b)] = xf j/N (x, Q 2 , b) - xf j/A (x, Q 2 , b) 
= 8vr A(A - l)3te ( ^"^H £ dxjp^ff i4 \^Q 2 ,x F ,t mhl ) 

/OO POO 
dz x / dz 2 p A {b,z 1 )p A (b,z 2 )e i ^-^ mN , (1) 

where fj/ A (x, Q 2 ), fj/^(x, Q 2 ) are nuclear and nucleon pdfs, ff^\f3, Q 2 ,xip, £ m in) 
are diffractive pdfs, 77 = ReA dl ^ / ImA dl ^ m 0.17 and p A (r) is the nuclear 
matter density. 

Numerical studies indicate that dominant contribution to the shadowing 
comes from the region of relatively (3 = Q 2 /(M 2 + Q 2 ) corresponding to 
rapidity intervals < 3 for which small x approximation used in the BFKL 
approaches is not applicable (these approaches lead also to the diffractive 
cross section with ap ~ 1.25 while the HERA experiments find a soft energy 
dependence corresponding to the effective Pomeron trajectory with otp ~ 
1.11). 

The uncertainties of the predictions are related to the shadowing effects 
resulting from the interaction of the hard probe with iV > 3 nucleons. Re- 
cently we improved the treatment of the multiple interactions, based on the 
concept of the color fluctuations and accounts for the presence of both point- 
like and hadron-like configurations in the virtual photon wave function [5]. 
This allowed us to reduce significantly the uncertainties: the difference be- 
tween two extreme scenarios (models I and II in Fig{T]) is < 20% for A ~ 200 
and much smaller for light nuclei. 

2. Where a non-linear regime sets in? 

In the leading log approximation one can derive a relation between the 
QCD evolution equations and the target rest frame picture of the interaction 
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Figure 1: Predictions of the theory of leading twist shadowing for quark and gluon pdfs. 
Difference in the predictions of two models of color fluctuation illustrates the range of 
uncertainties of the predictions. 



of small color dipoles with targets expressing it through the gluon density 
in the target, see [6] and refs. therein. Matching the behavior of the dipole 
cross section in the pQCD regime and of large size dipoles in the regime of 
soft interaction, it is possible to write interpolation formulae for the dipole - 
nucleon cross section for all dipole sizes and describe the total cross section 
of DIS at HERA. 

To determine how close is the interaction strength to the maximal allowed 
by the unitarity it is necessary to consider the amplitude of the dipole - 
nucleon interaction in the impact parameter space: 

r.M) = ^-J^y 2 J d 2 qe^A hN (s,t), (2) 

where A hN (s,t) is the elastic amplitude of dipole - nucleon scattering nor- 
malized to ImA hN = iShNCtotihN). The limit T h (s, b) = 1 corresponds to the 
regime of the complete absorption - black disk regime - the maximal strength 
allowed by the S-channel unitarity. 

The t dependence of the dipole - nucleon elastic scattering amplitude can 
be obtained from the studies of the exclusive vector meson production in the 
regime where QCD factorization theorem for the exclusive processes allows 
to express relate the t dependence of the amplitude to the t-dependence of 
the gluon GPDs. 

Combining this information with the information on the total cross sec- 
tion of the dipole - nucleon interaction allows us to determine Th(s,b) as 
function of the dipole size. A sample of the results for qq dipole -proton (nu- 
cleus) interaction which represent an update of the analysis of [7j is presented 
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in Fig|2j One can see from the figure that for small impact parameters T 
for scattering off a proton and a heavy nucleus are comparable provided one 
takes into account the leading twist shadowing for the gluons. This reflects 
an observation that gluon densities in nuclei and proton at b = are rather 
similar. However very few processes with proton are sensitive to b = rather 
than (b), while in nuclei T(b) is practically constant for a broad range of b. 
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Figure 2: Impact parameter distribution of qq dipole interaction with protons and lead 
nucleus. The models for the gluon shadowing are the same as in Fig(T] 



Since the probability of the inelastic interaction, Pi n {b) = \ — \\ — Y{b)\ , 
Pin > 3/4 already for Y{b) > 0.5. Also, in the case of the gluon induced 
processes Pi n {b) is enhanced by a factor 9/4 (provided it is sufficiently small). 
Hence the results presented in Fig|2] indicate that gluon induced interactions 
are close to the BDR for a much larger range of the dipole sizes. 



3. Post selection effect in BDR - effective energy losses 

It was argued in [8] that in the BD regime interactions with the target 
select configurations in the projectile wave function where the projectile's en- 
ergy is split between constituents much more efficiently than in the DGLAP 
regime. The simplest example is inclusive production of the leading hadrons 
in DIS for Q < 2p t (BDR). Interactions with the target are not suppressed 
up to p t ~ pt(BDR), leading to selection of configurations in 7* where longi- 
tudinal fractions carried by quark and antiquark are comparable. The energy 
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splits before the collision - post-selection. As a result to a first approximation 
the leading hadrons are produced in the independent fragmentation of q and 
q: 

D^\ z ) = 2 j\yD h q {z/y)^{l + (2y - l) 2 ), (3) 

leading to a strong suppression of the hadron production at xp > 0.3. 

In the case of a parton of a hadron projectile propagating through the nu- 
cleus near BDR effective energy losses were estimated in Ref . [9] . For quarks 
they are expected to be of the order of 10-r- 15 % in the regime of the onset 
of BDR and larger deep inside this regime. Also the effective energy losses 
are somewhat larger for gluons as the g — > gg splitting is more symmetric in 
the light cone fractions than the qg splitting. 

4. Leading hadron production in hadron - nucleus scattering 

Production of leading hadrons with p t ~ few GeV /c in hadron - nucleus 
scattering at high energies provides a sensitive test of the onset of the BDR 
dynamics. Indeed in this limit pQCD provides a good description the forward 
single inclusive pion production in pp scattering [11] . At the same time it was 
found to overestimate grossly the cross section of the pion production in dAu 
collisions in the same kinematics. The analysis of [O] has demonstrated that 
the dominant mechanism of the single pion production in the NN collisions 
in the kinematics which was studied at RHIC is scattering of leading quark 
of the nucleon off the gluon of the target with the median value of x for 
the gluon to be in the range x g ~ 0.01 -r- 0.03 depending on the rapidity 
of the pion. The nuclear gluon density for such x is known to be close to 
the incoherent sum of the gluon fields of the individual nucleons since the 
coherent length in the interaction is rather modest for such distances (cf. 
Figjl]). As a result the leading twist nuclear shadowing effects can explain 
only a very small fraction of the observed suppression [TJ] and one needs a 
novel dynamical mechanism to suppress generation of pions in such collisions. 
It was pointed out in [T2] that the energy fractional energy losses on the scale 
of 10 -T- 15% give a correct magnitude of suppression of the inclusive spectrum 
due to a steep fall of the cross section with xf which is consistent with the 
estimates within the post-selection mechanism. 

An important additional information comes from the correlations studies 
where correlation of the leading pion with the pion produced at the central 
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rapidities x [13] H3| which corresponds to the kinematics which receives the 
dominant contribution from the scattering off gluons with x g ~ 0.01 -j- 0.02. 
The rate of the correlations for pp scattering is consistent with pQCD expec- 
tations. An extensive analysis performed in [9] has demonstrated that the 
strengths of "hard forward pion" - "hard i] ~ pion" correlations in dAu 
and in pp scattering are similar. A rather small difference in the pedestal 
originates from the multiple soft collisions. Smallness of the increase of the 
soft pedestal as compared to pp collisions unambiguously demonstrates that 
the dominant source of the leading pions is the dAu scattering at large impact 
parameters. This conclusion is supported by the experimental observation 
[TO] that the associated multiplicity of soft hadrons in events with forward 
pion is a factor of two smaller than in the minimal bias dAu events. A factor 
of two reduction factor is consistent with the estimate of (9] based on the 
analysis of the soft component of 77 = production for the forward pion 
trigger. Overall these data indicate that (i) the dominant source of the for- 
ward pion production is 2 — » 2 mechanism, (ii) production is dominated by 
projectile scattering at large impact parameters, (iii) proportion of small x g 
contribution in the inclusive rate is approximately the same for pp and dAu 
collisions. 

A lack of additional suppression of the x g ~ 0.01 contribution to the 
double inclusive spectrum as compared to the suppression of the inclusive 
spectrum is explained in the post-selection mechanism as due to relatively 
small momentum of the gluon in the nucleus rest frame putting it far away 
from the BDR. 

It is difficult to reconcile enumerated features of the forward pion pro- 
duction data with the 2 — > 1 mechanism [TO] inspired by the color glass 
condensate model. In the scenario of [TO] incoherent 2 — )- 2 mechanism is ne- 
glected, a strong suppression of the recoil pion production is predicted. Also 
it leads to a dominance of the central impact parameters and hence a larger 
multiplicity for the central hadron production in the events with the forward 
pion trigger. The observed experimental pattern indicates the models [TO] 
which neglect contribution of the 2 — > 2 mechanism and consider only 2 — ^ 1 
processes strongly overestimates contribution of the 2 — > 1 mechanism to the 
inclusive cross section. 
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Figure 3: Distributions of the leading-power LO cross section for pp —> tt°-k°X at y/S = 
200 GeV in momentum fraction x a , for r)i = 2 (left) and r/i = 3 (right). We have chosen 
Pt,i — 2.5 GeV and integrated over 1.5 GeV < Pt,2 < Pti and various bins in 7y 2 . Units 
are arbitrary. The inserts show the corresponding normalized distributions on a linear 
scale. 

5. Production of two forward pions and double-parton mechanism 
in pp and dA scattering 

In Ref. P2] we suggested that in order to study the effects of small x gluon 
fields in the initial state one should study production of two leading pions 
in nucleon - nucleus collisions. Recently the data were taken on production 
of two forward pions in dAu The preliminary results of the studies of the 
reactions pp — > 7r°7r°+X, d—Au — > 7t°7t°+X, where one leading pion served as 
a trigger and the second leading pion had somewhat smaller longitudinal and 
transverse momenta [T?] US] . The data indicate a strong suppression of the 
back to back production of pions in the central dAu collisions. Also a large 
fraction of the double inclusive cross section is isotropic in the azimuthal angle 
A(p of the two pions. In order to understand the origin of the suppression and 
other features of the data we performed a study [19] which is summarized 
below. 

5.1. Forward dipion production in pp scattering 

it is instructive to start with the case of pp scattering. 

5.2. Kinematics of the 2-^-2 mechanism of two pion production 

We start with the leading twist mechanism 2 — > 2 in which a leading 
quark from the nucleon and a small x gluon scatter to produce two jets with 
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Figure 4: Comparison of the cross sections of LT (solid curves) and double parton mech- 
anisms (dashed curves). Inserts show the ratios of double parton and LT mechanisms. 

leading pions. In this kinematics x g < M 2 (im)/x q SNN- Production of two 
pions which together carry a large fraction of the nucleon momentum can 
occur only if x q is sufficiently close to one. The cross section is given by 
the standard expression for the 2 — y 2 processes. In Figj3] the results of the 
calculation of the integrand of the integral over x q for the LT mechanism are 
presented. They show that average value of x q for typical cuts of the RHIC 
experiments is very high. 

5.3. Double parton mechanism 

It is more likely that two rather than one quark in a nucleon carry to- 
gether x close to one. This suggests that in the discussed RHIC kinematics 
production the "double-scattering" contribution with two separate hard in- 
teractions in a single pp collision could become important. Hence though 
the discussed contribution is a "higher-twist", it is enhanced both by the 
probability of the relevant two quark configurations and the increase of the 
gluon density at small x which enters in the double-scattering in the second 
power. 

One can derive the expression for the double-scattering mechanism based 
on the analysis of the corresponding Feynman diagrams and express it through 
the new rank two generalized parton densities (GPDs) in the nucleons, see 
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[2U] and refs therein. The cross section can be written in the form 



— j — = [ dx a dx b dz c dx a >dx b >dz c > 

dp TA d Vl d PT , 2 d V2 TcRi nt J 



j2 ~afc— »c<i J2 'a'b'^rc'd' 

flh*))?W^)*W^W (4) 

Here f^{x a ,x a i) is the double parton distribution. If the partons are not 
correlated, it is equal to the product of the single parton distributions. For 
simplicity we neglected here correlations in the target as in our case x's for 
gluons are small. The dimensional factor ^R\ nt characterizes the transverse 
spread of the parton distributions in two nucleons as well as possible trans- 
verse parton - parton correlations and can be expressed through the integral 
of the product of two rank two GPDs (in many experimental papers and some 
theoretical papers this quantity is denoted as a e ff though it has little to do 
with the interaction cross section of the colliding hadrons). In the approxima- 
tion when partons are not correlated in the transverse plane one can express 
7rRi nt through the single parton GPDs, leading to 7rR^ nt as 34mb [21 J. This is 
larger than TrRf nt ~ 15 mb observed at the Tevatron, indicating presence of 
the correlations in the nucleon. In our numerical calculations we will use this 
experimental value. We find that for the RHIC kinematics the only trivial 
correlation due to the fixed number of the valence quarks is important while 
the correlation between x a and x a > remains a small correction if we follow the 
quark counting rules to estimate the x a > dependence off f a ,a'(%a, x a ') for fixed 
x a . The results of our calculation are presented in Fig. [4j They indicate that 
that the LT and double parton mechanisms are comparable for the kinemat- 
ics of the RHIC experiments. This provides a natural explanation for the 
presence of a large component in the pp — > 7r°7r° + X cross section measured 
in [TTl IT8] which does not depend on the azimuth angle between two pions, 
Aip. In fact the number of events in the pedestal comparable to the peak 
around A(p ~ tc which is dominated by the LT contribution indicating that 
the LT and double-parton contributions are comparable. 

Hence we conclude that the current experiments at RHIC have found a 
signal of double-parton interactions and that future experiments at RHIC 
will be able to obtain a unique information about double quark distributions 
in nucleons. It will be crucial for such studies to perform analyses for smaller 
bins in r\ and preferably switch to the analysis in bins of Feynman x. 
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(a) (b) (c) 

Figure 5: Three double partem mechanisms of dipion production. 

6. Production of two forward pions and double-parton mechanism 
in dAu scattering 

Let us extend now apply our results to the case of d-A scattering studied 
at RHIC. In this case there are three distinctive double-parton mechanisms 
depicted in Fig. |5j The first two are the same as in the pA scattering - 
scattering of two partons of the nucleon off two partons belonging to different 
nucleons (mechanism a), and off two partons belonging to the same nucleon 
of the target (mechanism b) [22]. The third mechanism, which is not present 
for pA scattering is scattering of one parton of proton and one parton of 
the neutron off two partons of the nucleus. Let us consider the ratio of the 
double-parton and leading twist contributions for dA and pp collisions 

r„ = r. + * + r.= ?^/™M (5) 

a LT {dA) (T LT (pp) 

The contribution to ta of the mechanisms (a), (c) is given by [22] : 

r c = T(b)a eff ;r a = 1, (6) 

where T(b) is the standard nuclear profile function (J d 2 bT(b) = A). Here 
we neglected nuclear gluon shadowing effect which is a small correction for 
the double-parton mechanism (cf. Ref.[12]) but maybe important for the LT 
mechanism where x g maybe as low as 10~ 3 due to the leading twist shadowing 
(see discussion below). 

For the central d-Au collisions Ta 2.2/m -2 and so r^jr c ~ 1/3 and 
Tb + fc ~ 4.4. The contribution (b) can be calculated in a model independent 
way since no parton correlations enter in this case. r c w r& for moderate 
rapidities where correlations between partons are not important. It reaches 
r c ~ 2rft/3 for the kinematics where only valence quarks contribute but 
x q + x q i 1. For the very forward region r c rj, since the kinematic 
constrain x q + x q > < 1 is not present in this case. As a result r^Au for small 
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Figure 6: Ratio of double-parton leading-twist contributions in dA —> 7r°7r°X. The inserts 
show the ratios r c /rb and r a /r c . 



b becomes of the order ten, see Fig. [6] (r^ changes from ~ 9 to ~ 12 for 
7rR 2 int = 15 4- 20 mb. 

Since the single inclusive pion spectrum for 772 ~ 2 4- 3 is suppressed by a 
factor of the order = 1/34-1/4 we find for the ratio of the pedestals 

in dAu and pp: 

Rpedestal = TaRa^) ~ 2.5 4- 4, (7) 

which should be compared with the experimental value of R pe destai ~ 3. 
Hence we naturally explain the magnitude of the enhancement of the pedestal 
in central dAu collision. 

If most of the pedestal in the kinematics studied at RHIC is due to the 
double-parton mechanism, the uncertainties in the estimate of the rates due 
to this mechanism and uncertainties in the strength of the suppression of the 
single inclusive forward pion spectrum at b ~ would make it very difficult 
to subtract this contribution with a precision necessary to find out whether 
all pedestal is due to double-parton mechanism or there is a room for a small 
contribution of the 2 — > 1 broadening mechanism as it was assumed in [23J . 

The suppression of the away peak originating from the LT contribution 
is due to two effects: (i) the gluon shadowing for x ~ 10 -3 and b < 3 fm 
and Q 2 ~ few GeV 2 reduces the cross section by a factor of about two (cf. 
Fig. 1), (ii) stronger effect of effective fractional energy losses due to larger 
x of the quark in the LT mechanism than in the double parton mechanism, 
leading to a suppression factor of the order two [19] . Combined these effects 
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result in a suppression of the order of four as compared to the single pion 
trigger, and overall suppression of the order of ten. This is pretty close to the 
maximal possible suppression which could be estimated as the probability of 
the "punch through" mechanism - contribution from the process where a 
quark scatters off one nucleon but does not encounter any extra nucleons at 
its impact parameter. Probability of such collisions at b ~ for interaction 
with Au nucleus is of the order 5 -j- 10% [23] . 

The preliminary data [17] are described well as a sum of the away peak 
of the strength < 1/4 of the strength observed in pp scattering, and the 
pedestal enhanced by a factor of three relative to the pp case, corresponding 
to reduction of the away peak relative to pedestal of the order of ten. 

7. Conclusions 

The LT twist nuclear shadowing can be estimated with a small uncer- 
tainty using correspondence with hard diffraction leading to the expectation 
of a large gluon shadowing at x < 10~ 3 . This influences significantly esti- 
mates of the proximity to the black disk regime, so that for small h they 
are reached for proton and nuclei at comparable energies. For b ~ black 
disk regime is effective for the RHIC energies leading to post selection effect 
(effective fractional energy losses) for the propagation of the leading partons. 
This is the only effect which currently allows to explain the suppression of 
the 2—7-2 mechanism for production of leading pions which dominates in 
the pp scattering including survival of forward - central correlations and sup- 
pression of forward - forward correlations. The processes of the production 
of two forward pions in pp, dAu scattering due to double -parton mechanism 
explain the pedestal observed at RHIC and make it difficult to search for the 
signal of 2 — > 1 processes. Magnitude of the away peak is consistent with 
the expected magnitude of the suppression of the 2 — > 2 mechanism due to 
LT gluon shadowing and fractional energy losses due to the post-selection 
mechanism. 

Also these experiments for the first time provide an opportunity to ob- 
serve correlations of leading quarks in protons. 
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